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Abstract 



We consider zero-temperature, stochastic Ising models a* with nearest-neighbor 
interactions in two and three dimensions. Using both symmetric and asymmetric initial 
configurations a^, we study the evolution of the system with time. We examine the 
issue of convergence of o"* and discuss the nature of the final state of the system. 
By determining a relation between the median number of spin flips per site u, the 
probability p that a spin in the initial spin configuration takes the value and 
lattice size L, we conclude that in two and three dimensions, the system converges to a 
frozen (but not necessarily uniform) state when p ^ 1/2. Results for p = 1/2 in three 
dimensions are consistent with the conjecture that the system does not evolve towards 
a fully frozen limiting state. Our simulations also uncover 'striped' and 'blinker' states 
first discussed by Spirin et al. [1], and their statistical properties are investigated. 

1 Introduction 

Consider the stochastic process a* corresponding to the zero-temperature limit of Glauber 
dynamics for a ferromagnetic Ising model with Hamiltonian 



Here (■) denotes a sum over nearest neighbor sites only and x, ?/ G Z'^, the li-dimensional 
hypercubic lattice. So cr* takes values in S = { — 1,+1}'^ , the space of infinite- volume spin 
configurations on Z*^. The process a* can be interpreted as modelling the nonequilibrium 
dynamical evolution of a classical Ising ferromagnet following a deep quench. 

Perhaps the most basic question that can be asked is whether the process ex* settles down 
to a limit as t ^ oo. Knowledge of this fundamental long-time property is a prerequisite for 
studying persistence properties [2, 3, 4, 5, 6, 7], and useful for the understanding of domain 
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formation and evolution, spatial and temporal scaling properties, and related questions (for 
a review, see Rcf. [8]), as well as aging phenomena (see, e.g., [9, 10, 11, 12, 13]). 

Given the fixed lattice type and dynamics under consideration, the convergence of cr* 
can depend only on dimension d and on the starting spin configuration cr". (Other lattice 
types were considered in [14, 15].) Rigorous results exist for symmetric a'^ (i.e., for all x, 
(j° = ±1 with equal probability, independently of all other spins, corresponding to a quench 
from infinite to zero temperature in zero external field) in d = 1,2; here every spin flips 
infinitely often [14, 15]. Numerical studies by Stauffcr [2] indicate that this result may hold 
up to (i = 4, beyond which a* could have a limit. There are few numerical studies [16] of the 
long-time behavior of cr* asymmetric 0"° (corresponding to a deep quench in nonzero external 
field) in dimensions greater than one. 

In this paper we present the results of numerical studies that examine the question of 
convergence of ci* for = 2 and 3 with both symmetric and asymmetric initial conditions. 
It is known that the magnetization per spin at zero temperature 

M(t)= hm Y: ^^(t)^ (2) 

\ xeAL 

where A^, is an L'^ cube centered at the origin, is independent of time in one dimension 
for any initial magnetization [17] . This time-invariance of the magnetization under Glauber 
dynamics at zero temperature does not hold in higher dimensions; indeed, it is easy to see, 
using the independence of the spins in a^, that in 2D at time t = the sign of dM/dt is that 
of the initial magnetization. For d = 2 our results indicate that the long time behavior of the 
symmetric case is unstable to small perturbations in the following sense: for any deviation 
from up-down symmetry in the initial state, cr* flows towards the uniform flnal state which 
is magnetized in the direction of the initial asymmetry. While this result is not surprising, it 
settles the question of whether the nonexistence of could persist for a nonzero range of 
asymmetry in . Our results for 2D are summarized in Fig. 1. For a 2D lattice of finite size, 
a non-uniform final state called the frozen stripe state is also possible. Such a configuration 
occurs when the spins are aligned parallel to the axes and cease to flip thereafter. 

In addition to the frozen stripe state, a finite size 3D lattice can evolve towards a 'partially 
frozen' state (called 'blinker states' in [16]) when a subset of the spins continue to fiip 
ceaselessly, without any change in energy, while the remainder of the spins do not fiip any 
more. Our results for the symmetric case in 3D agree with those of Stauffer [2]. For p = 1/2, 
we find, in agreement with [16], frozen stripe states and blinker configurations, to be discussed 
in the following sections. 

We also discuss the question of extrapolation of results on finite size lattices to those of 
infinite extent. 
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Magnetization Fiow at Zero Temeperature for 1 D and 2D ising iUlodeis 
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Figure 1: Sketch of magnetization flow at zero temperature under the Glauber dynamics 
discussed in the text, (a) One dimension, (b) Two dimensions. 



2 Definitions and Previous Results 

So far we have defined the problem informally: following a quench from infinite to zero 
temperature, in zero or nonzero external magnetic field, and subsequently evolving through 
standard Glauber dynamics in zero field, will the spin configuration eventually settle down 
to a final state, or will it continue to evolve forever (and if so, in what sense)? 

We now state the problem more precisely. The initial spin configuration (j° is chosen from 
a Bernoulli product measure (denoted P^'O) with each spin independently having value +1 
with probability p and —1 with probability 1 — p. The initial magnetization Af (0) = 2p — 1. 
The symmetric case corresponds to p = 1/2 and has M(0) = 0; any other P^o will be called 
asymmetric. 

The (single-spin-flip) dynamics will be taken to be continuous time, given by independent 

rate 1 Poisson processes at each x corresponding to those times tx (which can be thought of 
as clock rings at x) when a spin flip (o"*"'"'^ = — o"*"*^) is considered. If the resulting change in 
energy is negative, then the fiip occurs with probability 1, and if positive, with probability 
0. If there is a "tie", i.e., the resulting change in energy would be zero, then the fiip occurs 
with probabihty 1/2, as determined by a fair coin toss. We denote by the probability 
distribution on the realizations ou of the dynamics. 

So there are two sources of randomness, in the dynamics uj and in the initial spin config- 
uration (7°. The joint distribution of uj and (T° will be denoted P. All results to be discussed 
below are to be understood as occurring with probability one in P. 

The question posed in the introduction is whether cr* has a limit, with P-probability one, 
as i — > oo. This is equivalent to cr* flipping only flnitely many times for every x. Using the 
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same nomenclature as in [15], we call such an x an ^-site for finite); if cr* fiips infinitely 
often it is an X-sitc (X for infinite). By translation-ergodicity, the collection of JF-sites (resp., 
X-sites) has (with P-probability one) a well-defined non-random spatial density pjr (resp., 
px)- The densities pjr and pj {= 1 — pj^) can depend on both d and p. A particular model 
will be said to be of type JF, X or Al (for mixed) according to whether = 1, px = 1 or 
< PtiPx < 1, respectively. 

In one dimension the analysis is particularly simple, and it is fairly straightforward to 
show that, for any < p < 1, the model is type-X [14]. (The proof given there is for 
p = 1/2 only, but the result is not hard to extend to asymmetric models.) In [14] it was 
also shown that the homogeneous ferromagnet with p = 1/2 on is type-X; essentially 
the same argument holds for the homogeneous antiferromagnet on with symmetric initial 
conditions. No results were obtained for asymmetric initial conditions. 

We note for completeness that the ± J spin glass on with symmetric initial conditions 
was shown to be type Ai [19]; the proof is far more involved than those for the homogeneous 
cases. These are models in which couplings are chosen a priori from a symmetric Bernoulli 
distribution, with values +1 or —1, and are thereafter quenched. (In fact a much wider class 
of related models was also shown to be type-Al, but these are less relevant to the models 
studied here.) 

Under periodic boundary conditions, convergence to a final -|-1 or —1 state is not guar- 
anteed: the lattice can form domain walls parallel to the lattice axes. These walls are stable 
under the dynamics (hence referred to as the frozen 'stripe' states), so any state with one or 
more such walls (non- intersecting, i.e., all in either the x- or the ^/-direction) is a final state. 
A simple example of a domain wall in a 4 x 4 lattice is shown below: 

+ + + + 
+ + + + 



Investigations into this freezing phenomenon for a zero-temperature Ising ferromagnet 
having single-spin flip Glauber dynamics indicate that the system need not always evolve 
into the uniform state. Studies investigating the flnal state of zero-temperature Ising ferro- 
magnets [16] by Spirin et al. indicate that when d = 2, the system reaches either a frozen 
'striped' state (i.e., with one or more domain wall pairs) with probability ~ 1/3 or the uniform 
state with probability pa 2/3 [18]. Moreover their study indicates that these frozen states 
persist with positive probability as L ^ oo. When d > 2, persistent conflgurations other 
than the striped or uniform state are also observed. In these, the system wanders through a 
series of iso-energy partially frozen states. These metastable states are called blinkers in [16] 
and consist of localized sets of spins that flip forever with zero energy cost, surrounded by 
frozen spins. These blinkers become more numerous as the lattice size increases and the 
system is more likely to be trapped in one of these states. Presence of blinkers in the infinite 
system would indicate that it is of type-7W. These results will be discussed further in the 
Sec. 4. 
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Spin Flips as a Function of Initial Spin Probability for a 2D Ising lUlodel 
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Figure 2: In A/"/ InL, as a function of p in a 2i5 Ising Model. At least 10^ independent runs 
have been considered for L < 256, and at least 5 x 10^ independent runs for L > 256. 

3 Simulation Results 

Zero-temperature Glauber dynamics has been implemented by each spin independently cho- 
sen to flip through a Poisson process having rate parameter /i = 1. When a particular site 
is selected, it assumes the sign of the majority of its neighbors. In case of a tie, the sign of 
the spin is determined by a fair coin toss. 

Results depicting the median value of spin flips. A/", for 2D and 3D Ising models of 
different lattice sizes are shown in Figs. 2 and 3 respectively. For those runs in which 
the system converges to a final uniform state. A/" has been computed by averaging several 
independent (both in choice of initial spin configuration and in dynamical realization) runs. 

Our simulations in 3D also uncovered blinker states. When such a state appeared, the 
'unfrozen' spins (i.e. those that continued to flip) were tracked. 

From the simulation data for the 2D lattice, we observe that the number of spin flips 
is symmetric about a maximum at p = 1/2 irrespective of the lattice size. We also observe 
that the slope of this curve is sharpest near p = 1/2, while it remains relatively constant 
for very small or large values of p. For both the 2D and 3D models, the family of curves 
corresponding to different lattice sizes collapses to a single curve when In A/"/ ln2L is plotted 
vs. p. For very large lattice sizes, we have observed that the system almost always converges 
to a ±1 state even for small departures from p = 1/2. 

It is interesting to study the fraction of striped states as a function of p and lattice size 
L. This has been studied for a 2D model shown in Fig. 4. We observe that when p = 1/2, 
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Spin Flips as a Function of Initial Spin Probability for a 3D Ising IVIodel 
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Figure 3: In A/"/ InL, as a function of p in a 3D Ising Model. At least 10^ independent runs 
have been considered for all values of L, and only configurations that settle down to a fully 
frozen final state are included. 

as the lattice size increases, the fraction of the spins converging to a final striped state 
asymptotically tends towards 1/3, in agreement with the results in [16]. When p 1/2, the 
fraction of final spin configurations in the striped state decreases as lattice size increases; 
i.e., the uniform state is increasingly favored. No blinker configurations are present in the 
2D lattice. These results are consistent with the known result [15] that the infinite square 
lattice is type-X for p = 1/2, and they indicate that otherwise the lattice is type-jF. 

Fig. 3 shows results for the 3D cubic lattice; the data collapse similarly to the 2D case. 
As expected, the maximum number of spin flips occurs when p = 1/2. Fig. 5 shows the 
occurrence of striped states as a percentage of the total number of final states while Fig. 6 
depicts the occurrence of blinker states. Comparison of Figs. 5 and 6 clearly shows that, as L 
increases, the fraction of blinker states increases at the expense of the frozen striped states, 
in agreement with [1]. These trends suggest that for larger lattices, even small deviations 
from p = 1/2 causes the system to converge to the uniform state, while for p = 1/2, most 
final states are blinkers. In the latter case it seems likely that as L ^ oo, almost every final 
configuration is a blinker state (but sec the discussion in Section 4). 

The median number of spin flips per site, u = N'/L'^, was studied in 2D and 3D models 
for different system sizes and different p's (0.45, 0.50, 0.55), when the system converged to 
either a striped or a uniform state. For p = 0.45 and 0.55, u increases with L for small system 
sizes, but rapidly levels off after (roughly) L ~ 100 in 2D and L ~ 20 in 3D. Interestingly, 
these correspond to samples with approximately equal numbers of spins. Significantly, in 
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Occurrence of striped states in a 2D ising Model 
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Figure 4: Occurrence of striped states as a function of p for a 2D Ising Model. At least 10^ 
independent runs have been considered for p 7^ 1/2, and at least 10^ runs for p = 1/2. 
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Figure 5: Occurrence of striped states as a function of p for a 3D Ising Model. At least 10^ 
independent runs have been considered for p 1/2, and at least 10^ runs for p = 1/2. 
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Occurrence of Blinker configurations in a 3D Ising IVIodel 
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Figure 6: Occurrence of blinker states as a function of p for a 3D Ising Model. At least 10^ 
independent runs have been considered for p 7^ 1/2, and at least 10^ runs for p = 1/2. 



both 2D and 3D the rise of v with L continued steadily only at p = 1/2, where it increased 
monotonically as a linear function of L (cf. Figs. 7 and 8) up to the largest sizes considered. 
The rate of increase of u with L was, to within numerical accuracy, virtually identical in 
2D and 3D. This evidence indicates that, in terms of the overall classification of long-time 
dynamical behavior presented in Sec. 2, the Ising ferromagnet on behaves similarly for 
d — 2 and d — 3 as a, function of p. In particular, both are type-^ when p ^ 1/2, while 
for p — 1/2 the 3D lattice is either type-X as in 2D [14] or possibly type-^M. These data, 
however, do not reveal the nature of the final state(s) for large finite L; this will be discussed 
further below. 

We now study other aspects of the system evolution for p near 1/2. One important 
quantity of interest is the convergence time to a final state (for p 1/2) as a function 
of (p,d). Fig. 9 shows the result when p = 0.45, 0.50, and 0.55, for a 16 x 16 lattice. 
Times shown are in units of the inverse Poisson rate parameter (here set equal to one). 
Each simulation consists of 10^ independent runs where a final uniform or striped state was 
reached. 

For the 2D Ising model, when p = 0.45, 66.92% of the runs have final magnetization 
— 1. In a smaller percentage of cases (9.79%) the minority spin prevails, i.e., the final 
magnetization is +1. The remaining 23.29% are in the striped state. The average time to 
reach the —1 state is, not surprisingly, significantly smaller (37.64) than that taken to reach 
the +1 state (57.46). This extra time corresponds to large fiuctuations that take the system 
from the majority to the minority spin state. For the case when p — 1/2, the fraction of runs 
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Median value of Spin Flips per site (v) as a function of L in a 2D Ising Model 
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ure 7: Median value of spin flips per site, u, as a function of L in a 2D Ising Model. 



Median value of Spin Flips per site (v) and L in a 3D Ising Model 



p = 0.5 (Actual Data) 
p = 0.5 (Linear Fit) 
+ p = 0.45 (Actual Data) 
A p = 0.55 (Actual Data) 



-[•3.5 - 




ure 8: Median value of spin flips per site, u, as a function of L in a 3D Ising Model. 
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Histogram plots for a 16 x 16 lattice when p = 0.45, 0.50 and 0.55 
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Figure 9: Convergence time as computed by histogram plots for the 16 x 16 lattice when 
p = 0.45, 0.50 and 0.55. In all the plots the y-axis represents the frequency of occurrence. 
All the histograms are divided into 100 equally spaced bins and shown on a linear scale. 
The first row shows final magnetization states and the time taken to converge to the — 1 
and +1 states when p = 0.45. The second and third rows show the results for p = 0.50 and 
0.55 respectively. The results along the first and third rows are similar indicating symmetry 
around p = 0.50. 
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leading respectively to the uniform +1 and —1 final states are remarkably close (33.98% and 
33.76%, respectively) to the fraction 1/3. The system reaches a final striped state with a 
probability also close to 1/3, in good agreement with the results found by Spirin et al [16]. 
The corresponding time averages to reach these states are likewise very close (48.17 and 
47.67, respectively). The results for p — 0.55 mimic those of p = 0.45, with the relative roles 
of the +1 and —1 states switched. In this case 66.96% reach the +1 final state while 9.66% 
reach the —1 state. The corresponding time averages are 37.51 and 57.37, respectively. 

These tendencies are confirmed using larger lattices, as the next few figures show. In 
these, we compare results between a {2D) 64 x 64 lattice and a (3D) 16 x 16 x 16 lattice, 
which have the same total number of spins. The results are based on 10^ independent runs 
for either case. Graphs displaying final magnetization states and convergence times for the 
2 cases when p = 0.45 and p — 0.50 are displayed in Figs. 10 and 11 respectively. 

Comparing these two cases, we observe that for the 3-D Ising model, whenever convergence 
occurs, there is an overwhelming tendency for the majority spins to prevail in the final state. 
For p — 0.45, none of the 10^ runs in the 3D case converged towards the +1 state compared 
to 11 instances in the 2D case where the minority spin prevailed in the final state. For a 
64 X 64 lattice, 97.18% converged to the —1 state, 0.11% converged to the +1 state and the 
rest (2.21%) were striped states. Similarly, for the 16 x 16 x 16 lattice, 83.24% converged to 
the —1 state, 16.54% evolve into striped states and the remainder (0.11%) into blinker states. 
The average convergence time is smaller and more predictable for the 3D model than for 
the 2D case. For the 16 x 16 x 16 lattice with p — 0.45, the average time taken to converge 
to the —1 state is 32.00 compared to 304.53 for the 64 x 64 lattice for the same p. When 
p = 0.50, results are symmetric, as expected. Average times taken to reach the +1 and — 1 
states in the 16 x 16 x 16 lattice were 56.38 and 56.14, respectively, while they were 1729.43 
and 1837.02 for the 64 x 64 lattice. Compared to the 2D case, there is a smaller fraction of 
the spins converging to the final uniform state in the 3D case (2.9%) and this probability 
only decreases as L increases. Also when p = 0.50 in the 3D lattice, configurations that have 
not converged to the uniform +1 or —1 configurations tend to be closer to one of those two 
states than in the 2D lattice (cf. Fig. 10). 

Fig. 12 displays the nature of the blinker states in 3D. Typically, the number and position 
of blinker spins fiuctuates greatly from one configuration to another, hence the relatively large 
error bars. In these simulations, we considered at least 10*^ runs for L < 16 and 10"^ runs for 
L > 16. Two conclusions can be drawn from these data: first, at the lattice sizes studied the 
density of blinker spins is quite small; and second, there is a qualitative difference between 
the symmetric and asymmetric cases. In particular, when p = 1/2 the blinker density may be 
increasing with lattice size L, while it appears to remain fiat for asymmetric initial conditions. 
More work needs to be done for larger L to determine whether this trend persists. 
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Histogram plots comparing tlie final magnetization states for a 64 x 64 lattice and16 x 16 x 16 lattice 
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Figure 10: Histogram plots comparing the final magnetization states for a 64 x 64 lattice 
with a 16 X 16 X 16 lattice for p = 0.45, 0.50 and 0.55. In all plots the y-axis represents 
the frequency of occurrence while the x-axis represents the final magnetization states. All 
the histograms are divided into 100 equally spaced bins and shown on a linear scale. The 
symmetry for the p = 0.45 and 0.50 cases are evident from the first and last rows of the 
graph. 
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Histogram plots comparing time tal(en to reacli + 1 state for a 64 x 64 and 16 x 16 x 16 lattice 
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Figure 11: Histogram plots comparing the time taken to reach a final ±1 state for a 64 x 64 
lattice with a 16 x 16 x 16 lattice for p = 0.45, 0.50 and 0.55. In all plots the y-axis represents 
the frequency of occurrence while the x-axis represents the time taken. Note that the time 
bins are depicted on a logarithmic scale for the 2D lattice and linear scale for the 3D case. 
All the histograms are divided into 100 equally spaced bins. 
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Figure 12: Plots showing the bhnker density (that is, the fraction of bhnker spins in a given 
bhnker configuration) for a 3D lattice as a function of lattice size when p = 0.45, 0.50 and 
0.55. In these plots the y-axis represents the the mean number of blinker spins per blinker 
configuration, while the x-axis represents lattice size. 
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4 Discussion and Conclusions 



From our simulations for the evolution of the 2D and 3D Ising models with random initial 
conditions and evolving under zero-temperature Glauber dynamics, we draw the following 
conclusions: 

• For d = 2, the median number of spin flips is symmetric about a maximum at p = 1/2. 
The median number of spin flips per site, u, increases with lattice size but the rate 
of increase is significant only when p — 1/2 (cf. Fig. 7), when u — 0.12L + 0.33. At 

p — 1/2 and large L, the system evolves respectively towards the uniform +1 state, 
the uniform —1 state, and the striped state with equal probability 1/3 (over initial 
configurations and dynamical realizations). 

• However, it was proven in [14] that the infinite lattice with p = 1/2 is type-X. This 
may suggest that the L ^ oo limit, which distinguishes between striped and uniform 
states in a clear limiting sense, differs from the L = oo case. It is interesting to note 
that in the infinite square, at almost every sufficiently large time, any large region will 
locally resemble one of the two uniform phases or a domain wall phase. But this same 
region will continually approximate each of these three possible states infinitely often, 
separated (probably) by increasingly larger time intervals. On the infinite lattice, this 
state of affairs results from an unending supply of domain walls coming into any fixed 
region 'from infinity' — a state of affairs that cannot occur for a finite lattice, no 
matter how large. These results and considerations demonstrate that caution needs 
to be exercised in extrapolating numerical results to infinite lattices in two or more 
dimensions. 

• For d — 3, we again found the median number of spin flips to be symmetric about a 
maximum at p = 1/2. For this model the number of spin flips per site again increased 
linearly with L: z/ = O.IL + 0.33. It seems very likely that, as in the 2D case, a 3D 
lattice for p = 1/2 is not type-jF. The rate of increase of the median number of spin 
flips per site is, within the accuracy of the simulations, close to that of the 2D case. 

As L increases at p = 1/2, the fraction of final spin configurations in the blinker state 
increases rapidly (cf. Fig. 6). For the lattice sizes considered, however, the fraction of 
spins that are unfrozen in each of these 'blinker states' remains quite low. This is not 
surprising — it seems reasonable that every spin, or small connected cluster of spins, 
that cycles endlessly through a small number of states requires a significantly larger 
number of surrounding frozen spins to supply the 'boundary conditions' needed for this 
cycling to occur. 

It is unclear at this stage whether, as L increases, the fraction of 'unfrozen' spins in a 
typical final configuration levels off to a limiting value strictly less than (and probably 
considerably smaller than) one, as required for a blinker state, or continues to increase 
to one. The former case corresponds to type-A^ and the latter to type-X. Either way, 
these data are consistent with the conjecture that the symmetric case in 3D is not 
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type-^. Our data indicate that the symmetric case may be type-Al when L — oo; 
however, the previous example of the 2D case leaves open the possibility that the 
1/ — oo limit is tjpc-Ai while the L = oo case, i.e., the infinite cubic lattice, is type-X, 
and there is a corresponding dynamical singularity in the L ^ oo limit. 

• For p ^ 1/2, our results strongly suggest that the system is type-jT when L = oo and 
c? = 2 or 3, given that the average number of spin flips per site does not appreciably 
change with increasing lattice size. 

• Both L and d influence the probability of the minority spin prevailing in the flnal 
converged state. For any given p ^ 1/2, this probability decreases as L or d increases, 
thereby suggesting that the occurrence of a flnal minority spin ground state has zero 
probability when L — oo. 

• As might be expected, the average time taken to converge to a flnal majority state is 
smaller than that needed to converge to a minority state. For the same number of sites 
and any given p, dynamics on the 3D lattice shows a tendency to converge faster than 
in the 2D case. When p = 1/2, the times taken to reach the +1 or —1 flnal state is 
almost the same for both 2D and 3D. 
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